arXiv: 1508.0677lv2 [cond-mat.soft] 27 Apr 2016 


Percolation mechanism drives actin gels to the critically connected state 

Chiu Fan LecS 

Department of Bioengineering, Imperial College London, 

South Kensington Campus, London SW7 2AZ, United Kingdom 

Gunnar PruessneiHI 

Department of Mathematics, Imperial College London, 

South Kensington Campus, London SW7 2AZ, United Kingdom 
(Dated: April 28, 2016) 

Cell motility and tissue morphogenesis depend crucially on the dynamic remodelling of actomyosin 
networks. An actomyosin network consists of an actin polymer network connected by crosslinker 
proteins and motor protein myosins that generate internal stresses on the network. A recent dis¬ 
covery shows that for a range of experimental parameters, actomyosin networks contract to clusters 
with a power-law size distribution [Alvarado J. et al. (2013) Nature Physics 9 591). Here, we 
argue that actomyosin networks can exhibit robust critical signature without Hne-tuning because 
the dynamics of the system can be mapped onto a modified version of percolation with trapping 
(PT), which is known to show critical behaviour belonging to the static percolation universality 
class without the need of fine-tuning of a control parameter. We further employ our PT model to 
generate experimentally testable predictions. 

PACS numbers: 87.16.Ka, 64.60.ah, 05.65.-|-b, 64.60.De 


I. INTRODUCTION 


Actomyosin networks constitute an archetypal exam¬ 
ple of active matter mm- In a typical experimental 
setup, an actomyosin network consists of self-assembled 
actin filaments of variable length are connected via 
crosslinkers and myosins [3]. Actin filaments are two- 
stranded helical polymers of the protein actin, with a 
diameter of around 6 nm [4]. The protein fascin binds 
to two actin filaments and serves as a crosslinker to give 
connectivity to the whole network at high enough con¬ 
centration [2 . The molecular motor myosin II, generates 
internal stresses in the network by converting ATP to me¬ 
chanical forces [3]. In biological systems, myosins play a 
dominate role in muscle contraction and are responsible 
for the motility of eukaryotic cells. An actomyosin net¬ 
work can undergo contraction due to the internal stress 
generated by the myosins pulling two oppositely oriented 
actin filaments past each other HIS]. Such an active sys¬ 
tem has been shown to exhibit diverse patterns [am. Be¬ 
sides its intrinsic scientific interest, studying actomyosin 
networks is essential to our understanding of cell motil¬ 
ity, and cell and tissue morphogenesis mm- Recently, 
it was found that a planar actin network can be remod¬ 
elled by myosins in a way that exhibits robust scale in¬ 
variant structures similar to what is observed in static 
percolation at the critical point m- One of the most 
intriguing findings is that the characteristics of critical 
phenomena observed do not require fine-tuning of any 
control parameter, which suggests the presence of bona 
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fide self-organised criticality HIHISj- In addition to the 
experimental findings, the authors of [10| propose a de¬ 
tailed theoretical model that successfully accounts for 
the salient features of their experimental observations. 
Their model takes into consideration the various micro¬ 
scopic interactions among actin filaments, myosins and 
crosslinkers. 


In the present work, we asked a different question: 
what is the minimal physical model that explains the 
apparent self-organised critical nature of the active actin 
network observed? This question is strongly motivated 
by the expectations that at criticality, many microscopic 
details of the model can become irrelevant [Til ITS] . Be¬ 
low, we will argue that the underlying mechanism be¬ 
hind the observed critical behaviour of actomyosin net¬ 
work corresponds to a kind of percolation process known 
as percolation with trapping (PT) [T^, rather than self- 
organised criticality as envisaged by Bak, Tang, and 
Wiesenfeld [IIHS]. Our model not only recovers the 
most relevant experimental features, but allows in ad¬ 
dition to relate them to established results in the theory 
of percolation. 

The structure of the paper is as follows: In Sect.jnj we 
recapitulate the essential experimental findings in m- 
In Sect.jTTTjwe will argue that the underlying biophysical 
mechanism can be mapped onto a variant of percolation 
with trapping. In Sect.jTvjwe provide numerical evidence 
that the model indeed belongs to the static percolation 
universality class. We then derive predictions from our 
model, which include how the cluster size depends on the 
time of production, and the correlation between cluster 
sizes and their distances from the boundary. 
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FIG. 1: Schematics of the actin meshwork studied in the 
experiment by Alvarado et al. [10] (see also their Fig. 1(a)). 
Some actin filaments (black lines) connect to the top and 
to the bottom boundary. Interconnects are provided by 
crosslinking fascins (blue stars). Molecular motors myosins 
(red dumbbells) pull the filament together. The experiment 
was on a quasi-two dimensional plane of size 2.5 x 2.5mm^, 
and was observed over the course of about two hours. 

II. EXPERIMENTAL FINDINGS 

In their experiments, Alvarado et al. |I0j prepared a 
quasi-two dimensional network of actin filaments, about 
2.5 X 2.5mm^ in size, which are connected by fascin 
crosslinkers, illustrated schematically in Fig. Myosin 
motors pull the filaments together, thereby exerting 
contractile forces on the network [^. The resulting 
crosslinker unbinding and actin filament movement fun¬ 
damentally alter the underlying actin network (see sup¬ 
plementary movies in HO]). 

By tracing the trajectories of the actin filaments in the 
course of the contraction, the authors were able to recon¬ 
struct the connectivity of the actin network throughout 
the dynamical process. Specifically, depending on the 
number of crosslinkers in the system Me, three different 
regimes are observed m- 

1. At low Me, upon activation of the myosins, the 
actin network is contracted into many foci. By re¬ 
tracing the position of the actin molecules back to 
t = 0, one observes that the resulting foci origi¬ 
nated from actin clusters with areas of similar sizes. 

2. At intermediate Me, the retracing to the original 
areas indicates a power law distribution in sizes 
with an exponent of —1.91, which is similar to ex¬ 
pected exponent of —187/91 in random percolation 
at criticality m- 

3. At high Me, the actin network contracts to one 
single piece with the boundary of the experimental 
container discernible. 

The fact that the intermediate regime, in which scale in¬ 
variant size distribution is found, exists for a wide range 


FIG. 2: (a) The nodes adhered to the top and bottom bound¬ 
aries are red and with adhesion strength sa. The interior 
nodes (grey) have connection strengths drawn from a Pois¬ 
son distribution and the higher the strength, the darker the 
colour of the node. Each node is supposed to have a force / 
pulling on it on all sides and the weakest node will disappear 
first (b). 


of crosslinker concentration suggests that the critical sig¬ 
nature occurs without fine-tuning, as observed in systems 
displaying self-organised criticality. 


III. MODEL 

Alvarado et al. nni developed a model that incorpo¬ 
rates the microscopic details of the actin filaments un¬ 
der load, and the binding and unbinding kinetics of the 
crosslinkers. In the present work, we aim to account for 
the scale invariant nature of the experimental observa¬ 
tions by using a minimal model. The motivation behind 
this task comes from the general expectation that at crit¬ 
icality, many microscopic details become irrelevant. To 
obtain such a minimal description, we first list our key as¬ 
sumptions and our interpretations of the particular roles 
of the various ingredients in the experimental system. 


A. Initialisation of the actin network 

• We assume that the actin filaments form a planar 
network. For simplicity, we assume the network 
formed is a regular square lattice with N = L x L 
nodes. 

• The addition of crosslinkers serve to connect the fil¬ 
aments at the nodes in a random fashion. To model 
the randomness involved, we assign a strength value 
Sfe to node k that is chosen from a Poisson distribu¬ 
tion with mean s ^ Mc/N. This is to imitate the 
random number of crosslinkers connected at each 
node. 
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FIG. 3: (a) As random deletion of nodes start, the network 
starts to shrink, (b) As soon as a connected cluster detached 
from either top or bottom boundaries, it qnickly shrinks to a 
focal point or to the boundary, and as a result, node dele¬ 
tion is no longer possible (c). (d) Retracing the path of 

each node, the areas at the initial time corresponding to the 
three disjoint, fragmented clusters can be obtained, which are 
schematically shown in (d). 


• Contractile forces generated by the molecular mo¬ 
tor myosins induce pulling force fk on node k. Sim¬ 
ilar to the above, we assume that at the initial time, 
the pulling force is Poisson distributed with mean 
/ ^ Mra/N where is the number of myosin 
motors in the system. 

• Besides the pulling forces within the system, there 
is the adhesion that connects the actin network to 
the boundary. We assume here that the boundary 
adhesion strength is sa and following the modelling 
method in HO], we assume that the actin networks 
adhere to the top and bottom boundaries only. 
Note that although the adhesion strength here is 
in principle random, the exact values do not affect 
our analysis. 

The initial setup of the system is shown schematically in 
Fig-ia). We now describe the mysosin induced dynam¬ 
ics on the actin network. 


B. Dynamics 

• As the mysoin pulling force starts to act, we assume 
that the nodes rip (and thus effectively disappear) 
sequentially according to the reverse order of their 
unbinding propensity. The disappearance propen¬ 
sity of node k is given by a function B{sk, fk) such 
that B decreases with Sk but increases with fk- 


Physically, the disappearance of a node amounts to 
having the crosslinkers on that node unbind from 
the actin filaments at that node due to loading (Fig. 
Hb)). After a node disappears, the pulling forces 
are modified to achieve the force balance in the sys¬ 
tem [To]. This implicitly assumes that force bal¬ 
ance equilibration is fast compared to crosslinker 
unbinding uni. In other words, we assume that the 
unbinding times scale is much greater than the 
force balance equilibration time scale r/: ^ r/. 

• Besides the mysoin unbinding events, connected 
clusters of actin filaments can also contract, which 
we assume to occur whenever a connected cluster 
is detached from eit/ier boundary (Fig. |^. We also 
assume that the contraction time scale Tc is small 
compared to the unbinding timescale (t„ ^ Tc), 
so that once a cluster is free to contract, it contracts 
to a focal point almost immediately before the next 
unbinding event occurs. Once detached from the 
boundary the cluster freely contracts and all forces 
due to mechanical tension cease. The contracted 
cluster is thus removed from further dynamical evo¬ 
lution of the system, i.e., node deletions no longer 
occur in this cluster. 

Based on the above model ingredients, we can now ex¬ 
plain why the system is partitioned into three regimes as 
Me varies: if Me is small, many nodes will be empty and 
so the actin network would not be connected as whole, 
i.e., the network does not percolate. As a result, the 
actin network fragment into many small foci. On the 
other hand, if Me is high, then the average strength of a 
node may be higher than the adhesion strength (s > sa). 
As a result, the whole network detaches from the adhe¬ 
sion boundary due to the internal contraction, which is 
observed experimentally. The interesting regime is the 
intermediate range of Me, where the actin network perco¬ 
lates and the crosslinker unbinding occurs predominately 
in the interior. This is the regime we will focus on from 
now on. 

Instead of incorporating all of the above elements in 
the modelling, we will make the crucial simplification 
here by eliminating the dynamical elements in the sys¬ 
tem that are of time scale shorter than the time scale 
of node deletion in the network. In other words, we 
will ignore the force equilibration and cluster collapse 
steps. In addition, we assume that the pulling force ex¬ 
perienced by every node is identical at each time step. In 
this mean-field approximation, we are effectively assum¬ 
ing that although the pulling force at each node may vary 
and fluctuate as rupture occurs throughout the mesh, the 
spatial correlation of the force fluctuations remains short 
ranged so that upon coarse graining, a well-defined mean 
value exists. As long as the (initial) forces do not ex¬ 
ceed the adhesion strength sa to detach the mesh as a 
whole, and as long as they are strong enough to produce 
effectively random rupture throughout the system, the 
crosslinker unbinding events will produce random bond 
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deletion between cluster. The actual form of the unbind¬ 
ing propensity function B is therefore irrelevant as far 
as static properties are concerned. Indeed, the resulting 
fragmentation can be simulated as follows: 

Start with a fully connected squared lat¬ 
tice and delete interior nodes at random. 

If a connected cluster detaches from either 
boundary, then it is taken out the system and 
its configuration recorded. 

We will refer to this as the “rupturing model”. We note 
that a model based on random deletion of bonds in the 
context of active gels has also been proposed recently 
m- The main difference between our model and that in 
m is that our model incorporates the effect of boundary 
adhesion into the model description, which we believe 
is of experimental significance. In fact, we believe that 
both models are variants of the percolation with trapping 
model detailed below, where boundary conditions play a 
key role in the model dynamics. 

C. Relation to percolation with trapping 

The scheme above is a variant of (invasion) percola¬ 
tion with trapping |16] . Invasion percolation {without 
trapping) [TU] generates a single cluster by successively 
occupying more and more sites from an initial seed on an 
infinite lattice, thought to be occupied by the invader. 
Assigning initially a random number drawn uniformly 
and independently from the unit interval, the next site 
to be invaded is the one with the smallest such random 
number among non-invaded sites neighbouring invaded 
sites. Islands of non-invaded sites surrounded by invaders 
may appear in the process but disappear again as the 
invasion progresses. Invasion percolation with trapping 
suppresses these events, as no further invasion is allowed 
once an “island” of non-invaded sites is fully enclosed by 
invaders. As a motivation, consider extracting oil in a 
porous medium by pumping water (invader) in one spe¬ 
cific site m- The oil content in a particular pocket can 
be extracted as long as the pocket is connected to the 
extraction site. However, if the oil pocket is completely 
surrounded by water, the oil content will be locked in¬ 
side and thus be insusceptible to further displacement by 
water. 

Although originally motivated by invasion percolation, 
occupation by invaders does not need to be restricted to 
sites neighbouring to invaded sites. When occupation 
by invaders can take place spontaneously throughout the 
lattice, except inside isolated islands of unoccupied sites, 
this is known as percolation with trapping (PT) [TB] . 

The rupturing model as described above may be 
thought of as a version of PT in reverse. Nodes disappear 
in a random sequence, but no further nodes are deleted 
inside a cluster once it is detached from a boundary. The 
key difference between traditional PT and the ruptur¬ 
ing model above is this: While in traditional PT unoc- 
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FIG. 4: Cluster density distribution. Simulation of our PT 
model on a 200 x 200 square lattice produces a power-law dis¬ 
tribution of cluster density (n^) with an exponent close to 2 
(black curve). The red curve corresponds to the power law 
distribution with exponent —187/91 ~ —2.05 which corre¬ 
sponds to the exponent of static percolation on a 2D square 
lattice m- The inset figure shows the clusters generated from 
our model. 


cupied clusters stop shrinking as soon as they become 
detached from the boundaries, the rupturing (disappear¬ 
ance of occupied sites) stops as soon as an occupied clus¬ 
ter becomes detached from a boundary. In other words, 
our evolution algorithm is based on clusters of occupied 
sites becoming disconnected from the boundaries and we 
analyse the statistics of occupied sites, whereas the tra¬ 
ditional approach is based on unoccupied sites becoming 
disconnected from the boundaries, while still taking the 
statistics of occupied sites. 

A second difference, which numerically turns out to be 
insignihcant, is that in our model, clusters are already 
taken out of the system if they detach from any of the 
adhesion boundaries. Fig. This is to reflect the ex¬ 
perimental reality, as it was observed that clusters at¬ 
tached to only one adhesion boundary will contract to a 
the boundary quickly m- Although many big clusters 
therefore contract to the boundaries, smaller ones form 
earlier in the process as they get cut off from the bulk. 


IV. RESULTS AND DISCUSSION 

In this Section we will first show that our rupturing 
model reproduces the scale invariant nature of the experi¬ 
mental observations. We will then present several exper¬ 
imentally verihable predictions derived from our model 
simulations. 
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FIG. 5: Finite size scaling. Simulations of our PT model 
on different sizes of square lattice. The inset shows that 
cluster density distribution collapse onto a single curve as 
predicted by static percolation theory, which dictates that 
s'^Us = ^{s/L‘^-f ) where rl? is a scaling function and the expo¬ 
nents are: r = 187/91 and d/ =91/48 |17| . 


A. Static percolation universality class 


Numerical simulations of the original PT model have 
amply supported that PT belongs to the static percola¬ 
tion universality class m- In terms of cluster size distri¬ 
bution, Fig.j^shows that our modified PT also preserves 
the characteristics of static critical percolation on a 2D 
square lattice. As demonstrated in Fig. although the 
absolute value of the power law scaling obtained numer¬ 
ically seems to be smaller than the expected exponent of 
187/91 ~ 2.05, this is likely due to finite size effects. The 
collapse shown in the inset of Fig. [fusing exponents pre¬ 
dicted by static percolation theory clearly demonstrates 
that the modihed PT model presented here belongs to 
the universality class of static percolation m- We have 
also studied lattices with aspect ratios other than unity 
(the aspect ratio used in Figs. |4]|^ and found the same 
scaling behaviour, as expected from universality. 

We note that starting from a critically connected per¬ 
colation network, a no-enclaves version of the cluster size 
counting method has recently been proposed [IHlIin]- Us¬ 
ing this counting method, the power law exponent gov¬ 
erning the no-enclave cluster density as a function of clus¬ 
ter size is found to be lower than r ~ 2.05, which is thus 
closer to the experimentally measured exponent of 1.91 
m- In relation to our model, since the no-enclave mod¬ 
ification is ultimately based on a percolation network at 
criticality [21) . which is exactly the network generated 
by our model, the no-enclave modification can also be 
applied to our model straightforwardly. 


a) 
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FIG. 6: Predictions, a) Size of clusters as a function of the 
time of their free contractions. The inset shows an example of 
the temporal ordering (light to dark) of cluster contractions, 
b) Size of clusters as a function of the distance between their 
centres of mass and the closest adhesion boundary. Results 
are from simulations performed on a 200 x 200 square lattice. 


B. Predictions 

We have so far argued that for a range of experimental 
parameters, the contractile dynamics of an actomyosin 
network can be mapped onto a modified version of PT, 
which enables us to explain i) why critical behaviour is 
observed generically without fine tuning of a control pa¬ 
rameter, and ii) why the exponent observed in the clus¬ 
ter size distribution is close to that predicted by static 
percolation theory unj. We will now use our PT model 
to generate predictions that await experimental verifica¬ 
tions. 

Bigger clusters collapse later. Although the dynamical 
evolution of our model is purely sequential and physical 
time scales do not feature, we expect that the temporal 
trend of how the sizes of cluster produced vary as the sim¬ 
ulation progresses reflects what happens experimentally. 
As such, our first prediction is that within the critical 
domain (i.e, at the moderate range of crosslinker concen¬ 
tration Me), larger and larger clusters collapse onto foci 
as time progresses (Fig.j^a)). This prediction can be in¬ 
tuitively understood by the following argument: If large 
clusters collapse first, the network is unlikely to remain 
connected between the two adhesion boundaries, which 
would imply the end of the dynamical evolution since 
no nodes would be connected to both boundaries and no 
further rupture occurs. 

Cluster center of mass distribution. Fig. [^b) shows 
that on average, the centers of mass of larger clusters 
are further away from the nearest boundary, which can 
be trivially explained by the larger physical size of the 
clusters. A more interesting prediction here is that there 
is a peak in the curve shown Fig. [^b), which indicates 
that the centers of mass for large clusters stay away from 
the mid-line between the two adhesion boundaries. One 
plausible explanation is that large clusters are asymmet¬ 
ric in shape in such a way that a larger proportion of their 
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nodes are closer to the adhesive boundary that they re¬ 
main attached to. This asymmetry is built into the model 
by construction as a cluster is formed when it is discon¬ 
nected from one adhesion boundary while still being con¬ 
nected to the other. In other words, the cluster “rips at 
the thin end”, which implies that there is a thinned and 
a bulkier side. This predicted cluster shape asymmetry 
remains to be verihed experimentally. 

V. SUMMARY & OUTLOOK 

Contracting actomyosin networks can exhibit robust 
power-law distribution of cluster sizes [10]. The forma¬ 
tion of the clusters is due to the tearing apart of the 
initially homogeneous actin filament network by the con¬ 
tractile forces generated by the molecular motor myosins. 
Here, we provide a minimal model to account for the ro¬ 
bustness of the power-law scaling observed experimen¬ 
tally. Specifically, we argue that the key ingredients in 
the microscopic dynamics of an actomyosin network can 
be accounted for by a modified version of percolation with 
trapping. We show that this rupturing model exhibits 
critical behaviour that belongs to the static percolation 
universality class without fine-tuning. Furthermore, we 
generated specific predictions based on our PT model, 
which await experimental verifications. 

As far as the suggested connection to self-organised 
criticality is concerned, our model displays a form of self¬ 
organisation, as further rupturing stops when the clusters 
detach from the boundaries, i.e., close to the percola¬ 
tion threshold. This happens provided time scales are 
sufficiently separated. The scale invariant features are 
those found in static percolation. With a widely accepted 
definition of self-organised criticality in mind (such as 
non-trivial scale invariance, spatio-temporal power law 


correlation, self-tuning |13j l one might conclude that 
this really is a case of genuine self-organised criticality. 
However, there is no intermittency, no stationarity and 
strictly no dynamics. There is, in fact, no ongoing pro¬ 
cess that would enable self-tuning by some feedback loop, 
because scale invariance occurs as further evolution ter¬ 
minates. Our model therefore shares with invasion per¬ 
colation reservations that have been raised about it being 
an instance of self-organised criticality (35] [32 . 

In terms of future work, actomyosin networks consti¬ 
tute an archetype of active matter, which has been tradi¬ 
tionally modelled as a continuous medium with internal 
stresses generated within each volume element (33] [32 . 
The fact that the medium on which the stresses are trans¬ 
mitted can be transformed by the internal stresses, as in¬ 
vestigated here, challenges this traditional approach. In¬ 
deed, a theory of active matter on networks (rather than 
just lattices) would be required to properly account for 
the dynamical and mechanical properties of actomyosin 
networks in their full generalities. Another interesting 
future direction is to study the properties of myosins and 
crosslinkers as the network evolves. Here, we focused only 
on the cluster size distribution resulting from the active 
contraction and found that we recovered the static per¬ 
colation universality class. Active motility of particles in 
a continuous medium is known to lead to novel univer¬ 
sal behaviour (35] . The question therefore arises whether 
the dynamics of myosins and crosslinkers produces novel 
dynamical critical behaviour. 
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